NEW  YORK  UNtVERSlTY 

INarnVTE  of  MAH-fEMATlCAL  SCIENCES 
UBRARY 
AFCRC-TN-60-195  25  WaverJy  Place,  New  YoHc  3,  N.  Y. 

^^EET  PR/J 


■^^cccxx^'' 


"v^.         NEW  YORK   UNIVERSITY 

Institute  of  Mathematical  Sciences 
Division  of  Electromagnetic  Research 


RESEARCH     REPORT    No.    EM-156 


Diffraction  By  An  Elliptic  Cone 


LESTER  KRAUS  and  LEO  M.  LEVINE 


Contract     No.    AF     19(604)523 
MARCH,    1960 


^.-^^^-^•^'•. 


AFCRC-TN-60-195 


NEW  YORK  UNIVERSITY 

Institute  of  Mathematical  Sciences 

Division  of  Electromagnetic  Research 


Research  Report  No.  EM-I56 


DITFRACTION  BY  AN  ELLIPTIC  CONE 

Lester  Kraus 

and 
Leo  M.  Levlne 


^  yyi  X-*^ 


Leo  M.  Levine 


JM^h  ir. 


6,!^ 


Lester  Kraus 


M^  tt^  'Ke^'V 


Morris  Kline  Dr.  Werner  Gerbes 

Project  Director  Contract  Monitor 


The  research  reported  in  this  document  has  been  sponsored 
by  the  Electronics  Research  Directorate  of  the  Air  Force 
Cambridge  Research  Center,  Air  Research  and  Development 
Command,  under  Contract  No.  AF  19(60^)5258. 


Requests  for  additional  copies  by  agencies  of  the  Department  of  Defense, 

their  contractors,  and  other  Government  agencies  should  be  directed  to  the: 

ARMED  SERVICES  TECHNICAL  INEORMA.TION  AGENCY 
DOCUMENTS  SERVICE  CENTER 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 

Department  of  Defense  contractors  must  be  established  for  ASTIA  services 

or  have  their  'need-to-know'  certified  by  the  cognizant  military  agency 

of  their  project  or  contract. 

All  other  persons  and  organizations  should  apply  to  the: 

U.S.  DEPARTMENT  OF  COMMERCE 
OFFICE  OF  TECHNICAL  SERVICES 
WASHINGTON  25,  D.C. 


-  1 


Abstract 


The  boundary  value  problems  corresponding  to  a  spherical  wave  incident 
on  an  acoustically  soft  or  hard  elliptic  cone,  or  plane  angular  sector,  are 
solved.  The  method  of  separation  of  variables  in  a  sphero-conal  coordinate 
system  is  used.  The  special  functions  which  arise  are  simply  periodic  and 
nonperiodic  solutions  of  the  Lame''  differential  equation. 
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1.   Introduction. 

In  this  article  we  shall  derive  two  Green's  functions  for  the  reduced 
wave  equation  in  the  exterior  (or  interior)  of  an  elliptic  cone.   One  of 
them  and  the  normal  derivative  of  the  other  vanish  on  the  stirface  of  the 
elliptic  cone.   Both  satisfy  the  radiation  condition.   Thus  they  represent 
the  fields  produced  by  a  point  source  and  diffracted  by  a  soft  or  hard  cone 
respectively.   Since  a  plane  angular  sector  is  a  degenerate  elliptic  cone, 
the  problem  of  diffraction  by  a  plane  angxilar  sector  is  included  as  a  most 
Important  special  case.  These  solutions  can  be  used  to  determine  the  be- 
havior of  the  field  near  the  vertex,  this  being  carried  out  in  [ij .   They 
can  also  be  used  to  obtain  the  diffraction  coefficient  of  the  vertex  which 
is  needed  in  Keller's  geometrical  theory  of  diffraction  [2,3]- 

The  Green's  functions  are  constructed  by  the  method  of  separation  of 
variables  in  a  suitable  coordinate  system.  This  system.  Introduced  in  [^4-] , 
is  a  '\mif orroized '  version  of  the  sphero-conal  coordinate  system.   In  the 
analysis,  there  are  several  features  not  usually  encountered  in  separable 
problems.  First,  the  singular  siirfaces  of  the  coordinate  system  play  an 
important  role  in  determining  the  boimdary  conditions  for  the  resulting 
Stitrm-Llouville  problems.  Second,  the  two  separation  constants  occur  in  both 
angular  Sturm-Llouvllle  equations.  Therefore  they  must  be  determined  sim- 
iiltaneously .   In  this  respect  the  problem  resembles  that  of  determining  the 
vibrations  of  an  elliptical  membrane. 

The  special  functions  which  arise  are  simply  periodic  and  nonperlodlc 
solutions  of  the  I^me'' differential  equation.  The  doubly  periodic  Lame'' 
functions,  which  occur  in  classical  potential  problems  for  ellipsoids, have 
been  the  subject  of  numerous  investigations.   However,  the  simply  periodic 
Lame  ftmctions  have  not  been  studied  extensively  while  the  nonperlodlc  solutions, 
which  we  shall  call  nonperlodlc  Lame  functions,  have  apparently  not  been  studied 
at  all. 


We  will  formulate  our  diffraction  problems  as  follows.   Let  the  domain 
R  be  the  exterior  or  interior  of  a  semi-infinite  elliptic  cone,  or  the  exterior 
of  a  plane  angular  sector,  and  let  C  be  the  boundary  of  R.   We  wish  to  find  the 
two  Green's  functions  G^  "^(r^r")  and  G^^'(r,r')  defined  for  r,r '€R+C,  r  7^  r ' 
which  for  fixed  r'  satisfy  the  following  conditions  as  functions  of  r: 

(a)  ^^1^2)  ^  ^2^(1,2)  _^(^_^.) 

(b)  g(1^2)^,.(2)  .^^_^,^  g(1.2)^^(o)  i,j,^c-r' 

(c)  G^"*"^  =  0   on  C 

^(-t  (  2  ;  /  p  ^ 


bn 


0  and  grad  G^   exists  at  regular  points  of  C 


(1  2) 
(d)   G  '        satisfies  the  radiation  condition 


lim 
p  -*  0. 


/ 


^^    -   ikg(^'2)  ' 


dn 


centered  at  the  vertex  of  C  and  r  =  |r|.  The  fact  that  these  conditions 
determine  G  '    '   uniquely  has  been  proved  in  [5 J  . 

In  section  2,   we  describe  the  sphero-conal  system,  discuss  its  geometry, 
and  reformulate  our  diffraction  problem  in  terms  of  these  coordinates.   In 
section  5  we  apply  the  method  of  separation  of  variables,  discuss  the  resulting 
ordinary  differential  equations  and  derive  the  Green's  functions.  The  construction 
of  the  sphero-conal  coordinate  system  is  described  in  detail  in  Appendix  I,  and  a 
proof  of  the  completeness  and  orthogonality  of  the  Lame  products,  which  comprise 
the  angular  eigenfvmct ions  ,1s  given  in  Appendix  II.  The  proof  is  carried  out 
along  lines  indicated  in  \k\ , 

2.  The  reduced  wave  equation  in  sphero-conal  coordinates. 

Let  C  be  a  semi-infinite  elliptic  cone  with  vertex  at  the  origin  and  a 
directrix  consisting  of  an  ellipse  in  a  plane  x  =  const,  with  its  center  on 


-  2a 


source 


Figure  1 


The  domain  of  the  diffraction  problem  is  the  exterior  of  a 
semi-infinite  elliptic  cone.   The  elliptic  directrix  is  shown. 
■d     is  the  angle  between  the  generatrix  through  the  end  of  the 
minor  axis  and  the  positive  x-axis. 
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the  x-axls  and  its  major  axis  in  the  plane  y  =  0.   (See  Figure  1.)  Let 
R  be  that  domain  bounded  by  C  which  contains  the  positive  x-axis,  and  let 
4  be  the  angle  between  the  positive  x-axis  and  the  generatrix  through  the 
end  of  the  minor  axis  of  the  elliptic  directrix.  Thus  R  is  the  exterior  or 
interior  of  an  elliptic  cone,  depending  on  whether  o^'^  <n,  orO<'6  <- 
Let  6  be  the  semi-angle  subtended  at  the  origin  by  the  foci  of  the  ellipse  and 
let 

k  =  (1  +  tan^€  COST'S  Y^'^ 

(Here,  and  in  what  follows,  square  roots  are  always  to  be  taken  as  positive.) 
We  introduce  sphero-conal  coordinates  r,  •fl,  cp  by  the  equations 


2    2  1/2 
X  =  r  cos  ^(1  -  k'   cos  cp)  ' 


y 


r  sin  ■fl  sin  cp  (2.1) 


/,    .2   2.vl/2 
z  =  r  cos  cp(l  -  k  cos  ■0;  ' 

21/2 
Here,  k'  =  (l  -  k  )  '  ,  and  the  ranges  of  the  variables  are  given  by 

0  S  ^  s  rt 

0  ^  cp  ^  2n 

r  s  0 
The  construction  of  this  coordinate  system  is  described,  and  its 
orthogonality  proved  in  Appendix  I.   The  geometry  of  the  coordinate  surfaces 
may  be  described  briefly  as  follows  [see  Figure  2]  .  The  coordinate  r  is  the 
distance  to  the  origin,  so  that  the  set  of  coordinate  surfaces  r  =  const,  is, 
of  course,  the  set  of  spheres  centered  at  the  origin. 

The  coordinate  surface  -d   =  const,  is  a  semi- infinite  elliptic  cone 
whose  cross  section  in  a  plane  x  =  const,  is  an  ellipse  centered  on  the  x-axis, 
with  its  major  axis  in  the  plane  y  =  0.   ^  has  a  simple  geometrical  inter- 
pretation as  the  angle  between  the  generatrix  through  the  end  of  the  minor 
axis  of  any  cross  sectional  ellipse  and  the  positive  x-axis.  The  boundary 


Figure  2 

The  Sphero-conal  Coordinate  System.  The  sphero-conal  angular 
coordinates  ■&   ,   9  of  the  point  P  are  shown,  together  with  the 
coordinate  surface  ■3  =  -a,,  an  elliptic  cone,  and  the  surface 
cp  =  cp,  which  is  half  of  an  elliptic  cone.  The  cones  are  shown 
truncated  by  (cartesian)  coordinate  planes  through  P.  The  sphere 
though  P,  with  center  at  the  origin,  is  not  shown. 


cone  C  is  given  by  the  equation  ■&  =  -d    . 

The  surface  cp  =  const,  is  a  semi-infinite  elliptic  half -cone  whose 
cross  section  in  a  plane  z  =  const,  is  half  of  an  ellipse  centered  on  the 
z-axis^  with  its  major  axis  in  the  plane  y  =  0.  The  semi-ellipse  lies  above 
the  plane  y=OifO<cp<rt  and  below  the  plane  if  jt  <  cp  <  2jt.   The  coordinate 
cp  has  a  geometrical  interpretation  like  that  of  i3,  except  that  qp  takes  on  values 
greater  than  tt  when  the  generatrix  through  the  end  of  the  semi-minor  axis  (and 
the  semi-ellipse)  lies  below  the  plane  y  =  0. 

It  is  evident  that  the  coordinate  system  defined  by  (2.1)  becomes  the  familiar 
spherical  coordinate  system  when  k  =  1.  Thus,  when  the  semi-focal  angle  e  of  C 
vanishes, C  is  a  circular  cone.   Then  -d   becomes  the  polar  angle,  as  usual,  and 
the  semi-cones  cp  =  const,  become  half -planes. 

For  k  7^  1,  the  coordinate  surfaces  -a  =  0,  -d    =  n   are  degenerate  cones, 

namely,  plane  angvilar  sectors  in  the  plane  y  =  0.   Their  semi-angle  is 

arc  cos  k.   Similarly,  the  surfaces  corresponding  to  cf  =  0,  it,  2n  are  plane 

angular  sectors  in  the  plane  y  =  0  of  semi-angle  arc  cos  k ' .   Consequently, 

the  problem  of  diffraction  by  a  plane  angular  sector  corresponds  to  the  case 

^  =  Jt. 
o 

The  metric  coefficients  h  ,  h  ,  h  in  the  expression  for  the  line  element 

2     2    2     2    ?     2   2 
ds  =  h^  dr  +  h^  d'a  +  h^dcp 

are  easily  calculated  from  (2.1).  We  find 

2        2  2        2 

,2         2     k  sin  iS  +  k '   sin  qp  /•o  o^ 

1  -  k  cos  -d 


2        2  2        2 

,2  2     k   sin  ■a   +  k  '    sin  a 

h_  =  r 


5  ,        ,.,2„   2, 


ka 


Using  (2.2)  we  obtain  the  contravariant  form  of  the  gradient: 


/     2   2 
Vl  -  k  cos  •a 


A   n  |2   2 
vl  -  k'  cos  q 


^  =  ^r  ^  ^  ^. 


+  1 


/^ 


?   ?      ?   2 
r  Vk  sin  •a  +  k'  sin  Cf 


r  v/k 


(2.5) 


Here  1  ,  1„,  1  ,  are  imit  vectors  In  the  directions  in  which  the  indicated 
r'  ^'     cp' 

variables  are  Increasing. 

For  the  Laplacian  we  obtain 


Here 


r  r 


?   ?       2   2 
k  sin  ■d  +  k '  sin  q 


I  2   2~ 

Vl  -  k  cos  ^ 


hQ- 


■,2   2,   S 
1  -  k  cos  ■a  -r 


(s-*) 


(2.1f) 


yi 


^^(^ 


Before  formulating  our  diffraction  problem  in  sphero-conal  coordinates, 
we  must  derive  some  additional  'matching'  conditions,  which  arise  because  of 
the  peculiarities  of  the  coordinate  system.  The  first  of  these  is  the  familiar 
kind  of  condition  that  occurs  in  spherical  coordinates,  namely,  periodicity 
with  respect  to  the  azimuthal  angle.   It  is  evident  in  this  case  that  any 
function  u(r,'a,q))  which  is  a  single-valued  function  of  the  point  (x,y,z)  must 
be  periodic  in  qj  with  period  2jt,  if  we  extend  the  range  of  9  in  (2.1)  to 
-00  <  q)  <  00. 

We  note  next,  that  the  surfaces  *  =  0  and  ^  =  it  are  singular  coordinate 
surfaces,  since  every  point  therein  (except  on  the  edges  q)  =  jt)  is  described 
in  two  ways.  The  two  descriptions  are  (r,0,q))  or  (r,0,'2rt-q))  on  the  sector 
^  =  o  and  (r,rt,q))  or  (r,it,2rt-q))  on  the  sector  ^  =  n.   (See  Figures  3  and  ^•) 
Therefore  we  must  have  for  all  q),  r 


-  5a 


^=7r 


Figiire  k 

The  Intersections  of  the  ■d   and  cp  coordinate  surfaces  on  a 
sphere  aro\ind  the  origin  as  seen  looking  toward  the  origin 
from  the  positive  z-axis.  The  arrows  indicate  the  directions 
in  which  the  variables  are  increasing. 


u(r,0,cp)  =  u(r,0,23t-cp) 
By  taking  into  account  the  periodicity  with  respect  to  cp,  this  can  be  written 
as 

u(r,0,fp)  =  u(r,0,-cp)  (2.5) 

We  have  placed  oxir  'obstacle'  C  so  that  the  sector  -d  =  ic  Is   either  outside  the 
region  R,  or  forms  its  boundary,  if  C  is  a  plane  angular  sector.   Therefore  u 
need  not  be  the  same  on  both  sides  of  the  sector,  so  we  do  not  get  a  condition 
for  -d  =  Ti. 

We  also  require  that  grad  u  be  the  same  for  both  descriptions  of  the  sector 
^  =  0.   It  is  seen  in  Figure  h,    and  easily  verified  analytically  that  T  and  1 
reverse  direction  when  crossing  the  sector  •3=0.   Hence  we  must  have,  for  grad 
u  to  be  continuous,  (taking  into  accotmt  the  periodicity  in  cp), 


u^(r,0,rp)  =  -u^(r,0,-cp) 
u^(i-,0,cp)  =  -u^(r,0,-cp) 

The  second  equation  is  seen  to  be  a  consequence  of  (2.5)  and  hence  may  be 
omitted . 

There  are  also  matching  conditions  for  grad  u  at  the  lines  on  the  plane 
y  =  0,  common  to  adjacent  sectors.   Fro;n  Figure  h   we  see  that  on  opposite  sides 
of  the  line  •fi  =  0,  cp  =  0,  for  example,  the  unit  vector  1  has  perpendicxilar 

directions;  the  same  is  true  of  1  .   This  holds  also  for  the  other  three  such 

cp 

lines.  We  require  that  as  we  approach  the  line  •S  =  0,  cp  =  0,  for  example,  along 
either  of  the  sectors  •6  =  0  or  cp  =  0,  we  obtain  the  same  limit  for  grad  u.  For 
the  components  of  grad  u  perpendicular  to  the  plane  y  =  0,  we  get,  approaching 
the  line  ^  =  0,  cp  =  0,  along  the  sectors  cp  =  0  and  ^  =  0  respectively 

The  \init  vectors  1  ,  1„,  1  are  easily  computed  from  (2.1)  since  1„  ,  for 
r   u   cp  "0 

example,  has  the  direction  of 


(2.6) 


^  =  ■|^[lx(r,^,cp)  +ly(r,^,cp)  +  kz(r,-S,cpij 


For  the  components  parallel  to  y  =  0  and  tangent  to  the  sphere  r  =  const,  through 
the  point  we  get 

Taking  the  derivative  of  (2.5)  with  respect  to  cp,  and  using  (2.6)  and  the  con- 
dition of  periodicity  with  respect  to  cp,  we  get 

u^(r,0,0)  =  u^(r,0,0)  =  0 

The  first  two  of  the  last  three  conditions  can  now  be  rewritten  as 

^      u^(r,a,0)  -  u^(r,0,0)  ^     u^(r,0,cf )-u^(r,0,0) 


^  u^(r,0,O)  -  u  (r,0,0)  u  (r,0,t|))-ujr,0,0) 


Upon  evaluating  these  limits,  the  first  condition  turns  out  to  be  (k,k'  /  O) 
u^(r,0,0)  =  u^^(r,0,0) 

This  equation  is  trivially  satisfied  at  any  point  in  a  neighborhood  of  which 

both  u„,  u  exist,  and  at  which  one  of  the  mixed  derivatives  exists  and  is 

■d'      'p 

continuous.   The  second  condition,  however,  is  not  trivial:  we  have,  for 
k,k'  ^  0, 

k'2u^^(r,0,0)  +  k\^(r,0,0)  =  0  (2.?) 

In  essentially  the  same  way,  we  get  the  same  condition  on  the  line  ■6  =  0,  cp  =  rt. 


In  oiir  problem,  the  other  two  singular  lines  are  either  outside 
the  region  R,  i.e.  inside  the  'obstacle'  C,    or  form  its  edges.   In  the 
latter  case,  grad  u  is  not  necessarily  continuous  at  such  singularities 
of  the  boundary  surface;  in  fact  it  is  generally  infinite  at  an  edge. 
Thus,  in  either  case,  ve  do  not  get  any  matching  conditions  on  the  other 
two  lines. 

Condition  (2.7)  is  a  new  condition;  it  does  not  follow,  in  general, 
from  conditions  (2.5),  (2,6)  and  the  periodicity  in  cp.  For  example,  the 
fiinction  cos(^+cp)  satisfies  these  conditions,  but  does  not  satisfy  (2.7). 
However,  If  u(r,'6,cp)  is  a  solution  of  the  wave  equation  (or  potential,  or 
heat  equation)  then  (2.7)  is  automatically  satisfied.  To  see  this  we  write 
the  wave  equation,  using  (2.4)  for  the  Laplacian,  and  we  get 


1_  i_   ,  2  du^ 

2  ar  ^"^  b^'   ^  2,,  2 
r  r  (k 


^ ^—  (1-k^cos  •6)u  +(l-k'  cos  cp)u 

sin^^  +k'2sin2cp)  L  ^^  ^ 

+  /c^u  =  0 


2  2 

+  k  sin'dcos^u  +  k'  sincpcoscpu 

2   2      2   2 
Upon  cross-multiplying  by  k  sin  •a  +  k '  sin  cp,  and  letting  •a  -*  0,  cp  ^  0,  we 

see  that  the  expression  in  the  brackets  must  vanish,  and  we  get  (2,7).  Therefore 

we  need  take  into  account  only  the  matching  conditions  (2.5),  (2.6),  and  the 

periodicity  in  cp. 

We  can  now  formulate  our  diffraction  problem  in  sphero-conal  coordinates. 

We  seek  functions  G^"'"^(r,6,cp;  r',6',cp'),  G^  ^(r,^,cp;  r',i9',cp')  defined  for 

r,r'  g  0,  0  ^  -a,  ■d  '  S  ^  ,-«>  <  cp,  cp'  <  <»,  with  (r,-a,cp)  4   (r ',■3 ',cp '+2m0t), 

m  =  0,1,2,...,  satisfying  the  following  conditions:  [(r',^',cp')  fixed]: 

a)  The  wave  equation 


2^  i_  /^2  ^^^'^h   ^  k.     ^^q(1,2)_^2^(1,2)  ^  _  5(r-r')5(^-^')5((p-cp') 
r"  ^  V^  ^r     ;   ^2  ^2^ 

fl  2^ 

b)   Except  where  ■d   =  -d  ' ,    r  =  r',  and  cp  =  rp'+2mjt,  m  =  0,1, 2,  .  . ,,  G^  '  ' 

is  of  class  C^  ^  in  the  range  0^«<-a,r^O,  and  of  class  C^°' ,    0  $  ■&   ^ 


(2.8) 


c)  The  boimdary  conditions 
G^"'"^(r,^  ,Cp;r',i3Srp')  =  0 


(2.9) 


G^^^r,i3^,cp;r%^',cp')  =  0 


d)  The  radiation  condition 
o  2jt 


11m     r^      r       / 
r  -<■  00  ^        J^ 


hG 


(1,2) 


,(1,2) 


Sr 


(Jdcpd-a   =  0 


e)      The  matching  conditions 


G^^'^\r,^,cp;r' 


)   =  G^^'^^(r,'a,rp+2n;r',-a',cp') 


G^^^^^(r,0,cp;r',^',q.')    =   G^^' ^  V,  0, -rp;r ',^ ',p  • ) 


Gl^-'^^(r,0,(p;r%-S',cp')   =   -G^^' ^^  (r,  0, -cp;r  S^a  Sep  ' ) 
■a  cp 

In  the  above  A  is  given  In  (2.4')  and  a"  is  defined  by 

1   ,  ,  , 


/,  1,2   2,    /,  ,  ,2   2~    r^ 
Vl-k  cos  •6   yi-k'  cos  cp 


Then  r  cfdi^iiSdcp  is  the  element  of  volinne. 

5.      Green's   fiinctlong. 

We  consider  the  homogeneous   equation  corresponding  to   (2.8),    namely 


r  r 


(5.1) 


-  10  - 

2 
Setting  u=  R(r)v(i3,cp)  in  (j.l)  and  multiplying  by  —  ,  we  get 

As  in  the  case  of  spherical  coordinates,  we  write  the  separation  constant  as 
v(v+l): 

Then 

1^  (r^R*)  +  pV  -  v(v+l)l  R  =  0  "  (5.2) 

A  V  +  v(v+l)V  =  0  (5-5) 

The  solutions  of  (3.2)  are  the  spherical  Bessel  functions.   In  particular 
we  shall  employ  the  solutions 


2/^   ■^v+l/2('^^) 

(5.M 

^^-)  =  ^r   ^iil/2(-) 

The  first  function  satisfies  the  condition  of  continuity  at  the  vertex;  the 
second  satisfies  the  radiation  condition.  The  values  of  v  are  determined  by 
the  eigenvalues  v(v+l)  of  (3.5)  with  one  of  the  boundary  conditions  V('6^^cp)  =  0, 
or  V  (-3  ,cp)  =  0. 

Next,  equation  (3.5)  separates  as  follows.   Let  V(^,q))  =  e('a)  ^  (cp). 
Substituting  in  (3.3)  we  get  the  Sturm-Llouvllle  equations 

(l-k^cos^^)^/^  |^r(l-k2cos2^)^/^B'(^)l  +  [v(v+l)k2sin^-S-Hi]e(^)  =  0        (3.5) 
(l-k'^co8^cp)^/^^ri-k'2co82cp)l/2  ^'(q,)1  +  ^v(v+l)k•2slnV^I]  1  (cp)  =  0     (5.6) 


-  11 


Here  n  Is  the  separation  constant.  We  have  in  each  case  the  trigonometric 
fonn  of  the  Lame  differential  equation.  We  now  apply  our  matching  conditions 
to  the  eigenfunction  V(i3,cp)  =  F(^)  ^  (cp).   The  periodicity  condition  implies 
that  ^  (en)  must  have  period  2n.   Conditions  (2.11)  and  (2.12)  imply  that 

P(o)  $  (t)  =  e(o)  J  (-qp) 

9'(0)  J  M   =  -   P'(0)  I   (-cp) 

Equation  (5.?)  states  that  either  P(o)  =  0  or  ^  (cp)  =  ^  (-cp),  while  (5.8) 
states  that  either  P'(o)=  0  or  ^(rp)  =  -^(-cp).  We  cannot  have  both  R(o)  =  0 
and  '^'(O)  =  0,  since  this  woiild  mean  that  9{-d)   =   0.   Moreover,  since  ^  (cp) 
cannot  he  both  even  and  odd,  i.e.  identically  zero,  we  must  have  either 

e(0)  =  0     and     1  (<?>)  =  -  1  (-cp) 
or 

e'(0)  =  0    and     ^  (cp)  =  $  (-9) 

Since  (5-5)  ia  invariant  under  a  change  of  sign  of  the  independent  variable, 
we  see  that  if  P'(0)  =  0,  ^i-d)    is  an  even  solution  of  (5-5),  and  if  e(0)  =  0, 
then  P(^)  is  an  odd  solution.  Thus,  the  matching  conditions  imply  that  the 
factors  e(iS),  ^  (cp)  of  the  eigenfunction  V(^,(p)  are  either  both  odd  or  both 
even. 

We  now  define  solutions  ye(0;k,v,n),  yo(^;k,v,|i)  of  equation  (3-5) 
satisfying  the  initial  conditions 

ye(0;k,v,u)  =  1  yo(0;k,v,u)  =  0 

ye'(0;k,v,n)  =  0  yo'(0;k,v,^)  =  1 

We  shall  call  ye(■l3;k,v,^),  yo(i3;k,v,n)  the  even  and  odd  Lame  functions, 
respectively.   In  this  notation,  the  associated  solutions  of  (5.6)  are,  then, " 
ye(9;k'_,v,-n)  and  yo(cpik ',  v,-^)  . 


It  is  a  consequence  of  the  Stiinn-Liouville  theory  that  for  any  v  there 
is  a  discrete  set  of  values  p.  ,  Up,.-,   such  that  the  solutions  ye(cp;k ',v,-n^) 
of  (3.6)  are  periodic  with  period  2rt  (half  of  these  will  actually  have  period  n) ; 
similarly,  there  is  a  set  [x   ,    [i^,  . .  •   corresponding  to  odd  periodic  solutions 
with  period  2Tt.  Together,  these  sets  of  fimctions  comprise  a  complete  orthogonal 
set  over  the  Interval  [-«,«],  and  any  twice  continuously  differentlable  function 
can  be  expanded  in  a  uniformly  convergent  FoMxier   series '  of  these  functions 
(see  [9]  for  example).  These  functions  are  the  periodic  Lame  functions  and  have 
been  investigated  in  several  papers  (incei-  >(}.    gee  also  Erdelyi  ^ -I  ). 

On  the  other  hand,  the  functions of  •3  satisfying  one  or  the  other  of  the 
boundary  conditions  ye (•6  ;]s,  v,|i)  =  0  or  ye '(-3  ,k,v,\i)   =  0,  are  in  general  not 
periodic.   (The  case  •fl  =  n,  i.e.  when  C  is  a  plane  angular  sector,  is  an 
exception.  )  Again,  by  the  Sturm-Llouville  theory,  for  any  v,  there  is  a  sequence 
f  |i's  associated  with  even  solutions  forming  a  complete  orthogonal  set,  for 
nich  the  above  expansion  theorem  holds.  A  similar  statement  holds  with  respect 
to  odd  solutions. 

However,  we  have  here  the  two  parameter  eigenvalue  problem  of  determining 

all  the  pairs  (^q^'      f^^'^'      )   fo^"  which  the  solutions  ye(cp,k ',v^^'^%-|a^^' ^0 

d   2)   fl  2) 
are  2rt  periodic  and,  simultaneously,  the  solutions  ye(i3;k,v^  '      ,|j.^  '  ')   satisfy 

It  is  interesting  to  note  that  Ince  obtained  the  trigonometric  form  of 
the  Lame  equation  from  the  standard  elliptic  form  by  a  transformation  of  the 
independent  variable, for  the  purpose  of  simplifying  the  analysis;  on  the  other 
hand,  we  have  obtained  the  trigonometric  form  directly,  due  to  the  form  of  the 
coordinate  system  which  we  have  adopted. 

For  this  case  the  functions  •ye('6;k,v,^),  yo(^;k,v,n)  :are  periodic 
with  period  n,  2jt  or  4n.  A  detailed  discussion  is  given  in  [l]  . 
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the  boundary  conditions 

ye(^  jk,v^^^M^^))  =  0 
■'  ^  o'  '  en  '*^en  ' 

ye'(^  ;k,v(2^H^2^)  =  Q 
•'   ^  o'  '  en  '^en  ' 


(5.10) 


Similarly  we  seek  the  pairs  (v^  '   ,|i  '      )  for  which  the  odd  solutions  satisfy 


,(1.2)  Jl.2) 
on   'on 

the  same  conditions.   We  show  in  Appendix  II  that  the  Lame  products 

constitute  a  complete  orthogonal  set  of  eigenf unctions  of  (5.5)  for  each 
boundary  condition  . 

We  now  proceed  to  the  construction  of  the  required  Green's  functions. 
Set 

nal 


(5.11) 


(5.12) 


.1    .il'^hr)   yo(.;k,v(i'^),.(i'^))yo(.;k.  v(i'^),-.(i'^)) 
n=l 

We  substitute  (5-12)  into  (2,8)  and  Mse  (5.5)  to  simplify  the  resulting  equatloa.  Then  we 

multiply  by  c3ye(^;k,v^-^'^),^^-'"'^^  )ye(cp;k ',v^''"'^S-|i^-'"'^^),  and  integrate  over 

the  region  -n  S  q)  S  tx,    0  ^  ■&   ^  ■&    .      Using  the  orthogonality  of  the  Lame  products 


with  respect  to  the  weight  function  o,  we  get 


dR 


^1,2) 


1_  d^  (2  

2  dr  V^    dr 


^(1,2)^^(1,2)  ^^)- 


,(1,2) 


(5.15) 


ye(.';k,v(^>^),^(^-^))ye(cp';k^v(^-^),-^(^-^))6(r-r^) 


Ik 


-Jt  o   - 

(l   2) 
e  function  R\  '    '    satisfies  the  spherical  Bessel  eqxiation 

-j  +  rA2_^(l,2)   (1,2)^   1^(1,2)^^ 
/    [_       em   ^  em       'J  em 


cJd-adcp 


dr  V'''     dr 


except  at  r  =  r'.  There  it  is  continuous,  with  a  jump  in  the  derivative 
given  by 


(3.1M 


(3.15) 


dR 


(1,2) 


dr 


ye(.';k,v^l'2)   (1,2))   (  ^     (1,2)    (1,2), 
^    ^      ^  ^  em   ^"^em   i^    ^^    >      ^  em   '      em 


(3.16) 


.(1,2) 


Since  R^^  Mr)  is  to  be  continuous  at  the  origin,  it  must  be  proportional  there 
em 

to  j  /    Akt) .     To  satisfy  the  radiation  condition  it  must  be  proportional  to 

em 
h  -^  p'\('^^)  ^°^  large  r.   From  the  continuity  at  r',  ve  get 

em 

«^^'"  =  =i^''\(l,2)(-<)\(l,2)(->> 

em         em 

Here  we  have  used  the  common  notation 
r  =  min(r,r ' ) 

r  =  max(r,r  ' ) 
When  we  apply  the  jump  condition,  (5.I6),  to  the  left  side  of  (5.I7),  we  get 

ye(.-k,v(l^2)   (1,2)  (1,2)    (1,2) 


(3.17) 


(1,2) 


em     em 


(3.18) 


/<N 


(1,2)^,2 


^(J  (1,2)'  "^   (1,2)^ 
^   em       em    - 

Here  W(j,h)  is  the  Wronsklan.   We  have 


Kr ' 


W(j 


(1,2)'  \(1,2)^ 


(5.19) 


{Kr^Y 


-  15  - 

(l  2) 
We  proceed  in  exactly  the  same  way  to  find  the  coefficients  R^  '    ' {v)   of 

the  odd  Lame  products. 

Substituting  for  Rg^'^^Cr)  and  Rq^'^^(i-)  in  (5.12),  we  get  for  the 

Green's  functions (p=l, 2) 

00 

n=l  en    en        en 
ye(.A,v(P),.(P)),e(.M.,v(P)„(^))ye(,ik.,v(P),-.(P))ye(..;..,v(P',-,(P)) 

n=l  on    on        on 

Here,  N^"^'^^  is  defined  hy  {^.ih),   and  r^^'^^  is  defined  similarly  by 
en  on 

-li'^'  =  /  /  h^^^-(i.^),.(i'^))yo(.;..,v(i.^),-.(;'=))]'  .dM.       (,..1, 

-n  o 

The  expression  (3.20)  for  the  Green's  functions  is  our  final  result. 
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Appendix  I.  The  sphero-conal  coordinate  system. 

By  analogy  to  confocal  (ellipsoidal)  coordinates,  we  consider  the 

set  of  'confocal  cones' 

2     2      2 
f-  +^  +^—=   0  0  <  b^  <  c^  (l-l) 

5    l-b^   l-c^ 

These  surfaces  are  two-napped  elliptical  cones  with  a  common  vertex  at  the 
origin.  Exactly  as  in  the  case  of  confocal  coordinates,  we  can  easily  prove 
that  if  two  distinct  surfaces  of  (l-l)  have  in  common  a  point  other  than 
the  origin,  they  are  orthogonal  at  that  point.  To  do  so  we  define  F  by 

2     2      2 
F(x,y,z;|)  =  ^-  +  "^  +  -^ 

Let  (x  ,y  , z  )  be  a  point  common  to  the  surfaces 

F^(x,y,z)  =  F(x,y,z;|^)  =  0 

F2(x,y,z)  =  F(x,y,z;|2)  =  0  ^1  '^  ^2 

Then  we  have,  at  (x  ,y  , z  ), 

hF^  hF^       bF^   aFg   BF^  BFg      1,    r  -i 


This  proves  the  orthogonality. 

If  x,y,z  are  fixed,  (l-l)  becomes  a  quadratic  eq\zation  in  |  and  con- 
sequently has  at  most  two  distinct  roots.  Thinking  of  F(x,y, z;|)  as  a  fvmction 

2  2 

of  5,  we  easily  see  that  it  has  at  least  one  zero  v  between  0  and  b  ,  and  a 

2  2      2 

zero  |i  between  b  and  c  .  Therefore  it  has  exactly  one  such  zero  in  each  of 

these  intervals.   Hence  through  any  point  (x,y, z)  ^   (O, 0,0)  there  are  exactly 

two  cones 


-  17' 


2      2  2 

^  +—2: 5—  =  0 

2  ^  2  ^2  2  2   ^ 

H     n  -  D  C  -  |1 

2      2  2 

~2  "  -  2  2  2  2  "  ° 

V    b  -V  c  -u 


We  shall  use  these  two  sets  of  cones  as  coordinate  surfaces.  As  our 
third  set  of  coordinate  surfaces,  orthogonal  to  the  two  sets  of  cones,  we 
take  the  system  of  concentric  spheres  centered  at  the  origin.  Thus  we 
consider  the  transfonnation  of  coordinates  from  (x,y,z)  to  (r,|i,v) 

n2222  v.1,^^ 

a)r=x+y+z  c>b>0 

2  2  2 

b)      ^  +  -|— 2  -  -^  =   0  0  g  v^  S  b^  i  H^  S   c^ 

[i  |i   -b  c   -n 

2  2  2 


.)    ^ z. 

■I  o  o 


2       ^2      2  2     2^ 

V  b   -V  c   -V 

2  2  2 
These  equations  give  \l   ,v   ,t     in  terms  of  x,y,  z,  as  solution  of  quadratic 

equations. 

2      2 
To  get  the  inverse  transformation,  we  use  the  fact  that  \i     and  v  are 

the  roots  of  the  equation  ( in  | ) 

2     2      2         f. 

2L-  +-L_  +_L_^  Eli] 0 

i  l-b^   i-c^   5(5-b2)(5-c2) 

We  note  that  v{i)   is  a  quadratic  expression  in  I   with  leading  coefficient 

2    2    2    2  2   2 

X  +y  +z  =r,  and  roots  \i   ,    v    .      Hence 


^    l-b^   i-c^     5(|-b2)(^-c2) 
Upon  mioltiplying  (l-5)  by  |,|-b  ,5-c  ,  and  setting  |=0,  |i»b  ,  5=c  , 
respectively,  in  the  resulting  identities,  we  get  the  inverse  transformation 
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2  2  2 
r  V  |a 


^)    2  ^  r^(K"-b^)(b^-v^)  2    2    2  d-'^) 


,  2,  2  ,  2v  t  g  n   S  c 

b  (c  -b  )  "^ 


2/2  2w  2  2. 
2  _  r  (c  -[X    ){c   -V  ) 


2/  2  ,  2v 
c  (c  -b  ) 


For  a  given  set  of  values  r^n,v  the  siorfaces  (l-2),  i.e.  the  sphere  a) 
and  the  two-napped  cones  b)  and  c).  Intersect  in  eight  points,  in  general. 
We  attempt  to  achieve  a  1  -  1  correspondence  between  the  points  (x,y, z)  and 
the  coordinates  (r,|a,v)  by  proceeding  as  follows.  First  we  restrict  (r,|a,v) 
to  the  ranges  rSO,  b^^i^c  and  -b  S  v  S  b.  Then  r,\i,v   are  iiniquely 
determined  by  x,j,z.     To  determine  x^y, z  from  r,n,v  we  must  specify  which 
square  roots  are  to  be  taken  in  solving  (l-^).   In  (l-'^)  we  take  that 
sq\iare  root  which  has  the  same  sign  as  v . 

We  cannot  remove  the  remaining  ambiguities  in  the  signs  of  y  and  z  in 

/~2     2 
the  same  way.   Therefore  we  introduce  instead  of  ^i  the  new  variable  X  =  \/c   -^    . 

/~2  2        /~2     2 
We  assign  to  X   the  range  -  \/c   -h     ^  X  S  l/c  -b  .   Now  we  define  z  by  taking 

that  square  root  in  (l-^c)  which  has  the  same  sign  as  X. 

In  order  to  specify  which  root  is  to  be  taken  in  (l-^b)  it  is  convenient 

to  replace  v  andX  by  the  two  new  variables  •6  and  cp,  defined  by 

V  =  bcos'3  0  ^  ^   S  rt 

X  =    \/c^-ii       =     /c^-b^       cos  cp         0  ^  cp  S  2jt 


/  2     2  /  2     2 

\/[i  -b       =     \/c  -b        sin  (p 

Now,  upon  setting  k  =  -,  we  get  from  (l-h)   the  formulas  (2.1).   These  formulas 
give  the  desired  1-1  correspcnBience  b«tw»«ii  the  coordinates  (x,y,z)  and  (r,i8,qp) 
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except  for  the  singular  coordinate  surfaces  -d  =  0,   -d  =  it   (which  are  discussed 
in  section  2) . 

Upon  substituting  for  v  and  -  in  (l-2c)  we  see  that  the  surface  -d  =  ■& 
is  a  single  nappe  of  the  cone 


2             2 
_Z +  5 

sin^^^   (^  -  l)+sin^'fl^ 


The  nappe  of  the  cone  encloses  the  positive  or  negative  x-axls  depending  on 


/ 

+ 

2 
z 

.lt.n\ 

2        2                 ""l 

t^- 

■  1) 

^           °       cos^^ 
o 

If  (1-5) is  the  boundary  cone  C,  and  x  =  x  the  plane  of  its  directrix,  the 
latter  is  given  by  (I-6).   If  e  is  the  semi-focal  angle  of  C,    then  x  tane 
is  the  semi-focal  distance  of  the  directrix,  and  from  (I-6)  we  get 


>K      V  v2 


This  gives  us  k  in  terms  of  the  boundary  cone  C. 


Appendix  II.   Orthogonality  and  completeness  of  the  Lame'' products . 

Let  V  (■a,cp),  V  (■6,cp)  be  any  two  of  the  Lame'' products  of  (3-11),  satisfying 
the  same  boundary  condition.  We  shall  first  prove  that,  unless  V^  h  V^/  they 
are  orthogonal,  i.e. 
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(ii-i; 


First,  we  note  that  if  one  of  the  products  V  ,  V  is  odd  and  the  other  even, 
(ll-l)  follows  immediately;  for  then  V^VpO'woiLLd  be  an  odd  fimction  of  9,  and, 
integrating  first  with  respect  to  cp,  we  have  the  result.  We  next  assume  that 

V,  =  ye(^;k,v^^^^^^))ye(cp;k',v^^^-^^^^) 

1  J    \    '    '    em  '  em.^    ^    '      '    em  -^   em 

V^  =  ye(^;k,v^^^H^^^)ye(9;k',v(^V^i^^^) 

2  ^    ^    '    >    emg'^em^'-^  ^^'  em^'   "^em^' 

(The  proof  proceeds  in  exactly  the  same  way  if  V,,"^^  are  odd  Lame'' products. 

From  (5.5)  we  get  (lI-2) 

d^ 


or  if  they  satisfy  the  second  boundary  condition.)  Let  X.    =  v^    ' (v^        +  l); 
"^       ''  i    em   em 


i,  (d  -  A.sW^-[e/,  -  e,e;j|.  [(.^-x^^^s.n^.  .  .^  -  .(^)]— i 


1-k^cos^^ 


Upon  integrating  (lI-2)  from  0  to  6  ,  we  find  that  the  left  side  vanishes, 

since  the  quantity  inside  the  square  bracket  vanishes  at  both  end  points.   This 

is  so  at  i8  =  0  because  the  functions  are  even  and  at  ■d  =  i8  because  we  are 

o 

assuming  the  first  boimdary  condition.   It  is  also  clear  that  the  square  bracket 
vanishes  at  the  end  points  for  the  various  other  possibilities,  i.e.  for  a 
rigid  boiindary  or  odd  ftmctions  or  both.   Consequently  we  get 

(X-X,)k2  /    ^^       d^  =  (u^^^   -u^^b      /     .   ^  II-3) 

^12^    J  j ^ ^      ^^em    ^em^  J  1 ^ 2" 

o     \/l-k  cos"^^         "^  ^     o  Vl-k^cos^'fl 

Similarly,  upon  using  (5.6)  we  get 
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-n   \/l-k'  COB  cp  -jt   Vl-k'  cos  cp 


Now,  multiplying  opposite  sides  of  (lI-5)  and  (ll-^),  and  combining  terms 
we  get 


•It  o 

If  neither  n^"^^  =  ^l^"^^  nor  X,  =  X, 
"^em^    em  2       1    2 

(Il-l).   If  [1^-^'   =  M^-^^  but  X-  4  ^n>    then  the  integrals  on  the  left  sides  of 
em.^    em       1  '   2' 

(II-5),  (11-^)  vanish,  and  hence  the  double  integral  in  (II-5)  vanishes. 
Similarly,  if  X     =  X  ,  but  \i^        7^  ^   ,    then  the  integrals  on  the  right  sides 
of  (II-3),  (11-^)  vanish,  and  again  the  double  Integral  in  (II-5)  vanishes. 
This  completes  the  proof. 

We  show  next  that  the  lame'  products  (5 .11)  form  a  complete  set  of  eigen- 
fxinctions  for  equation  (JO)-  First,  we  recall  that  there  exist  complete '  sets 
of  eigenfunctions 

n  =  1,2,...  for  the  first  and  second  boundary  value  problems  respectively. 
This  is  proved  in  [5] .  We  shall  show  that  every  eigenfimction  V^  '    '   can  be 
expressed  as  a  finite  simi  of  Lame'  products,  and  hence  that  the  latter  constitute 
a  complete  set. 

Again,  for  def inlteness,  let  us  consider  a  particular  boundary  value 
problem  say  the  problem  of  the  rigid  boundary,  i.e.  ^/dn  =  0  on  C.   We 
shall  omit  the  superscripts  indicating  this.  First,  we  note  that  if  V  (iS,q)) 
satisfies  (5.5),  then  so  does  V  (iS,-cp).  But 


uj.on  \]j.^).      rirtju,  we  rectij-X  T^nttt  ^iisre    ex±tst   compxeue  seu! 
I  n  r  '      |v^^7  ^ith  associated  eigenvalues  1^^  (l^n  I' 
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Since  the  first  bracketed  term  is  a  solution  of  (3-5)  which  is  even  in  t),  and 

the  second  a  solution  odd  in  cp,  we  can  assume  without  loss  of  generality  that 

V  is  either  even  or  odd  in  cp. 
m 

Let  us  assimie  that  V  (■0,<p)  is  even  cp.   Let  v  be  its  eigenvalue.  Then  by 
the  remarks  following  (3-9),  V  can  be  expanded  in  an  infinite  series  of  periodic 
Lame'  fimctions  ye(cp;k,v  ,-\x.),    J  =  1,2,...   .  Because  of  the  evenness  of  V  ,  the 
Fourier  coefficients  of  the  odd  periodic  fimctions  vanish.   Let  us  normalize  the 
functions  ye(cp;k,v  ,-n .),  denoting  the  corresponding  normalized  fimction  by 
^  (cp).   Thus,  we  have 


/ 


'-^ ^^  =  ^ij 


/    2   2 
Vl-k '  cos  cp 

Consequently  we  have 

V^(^,cp)   =      ^    A.(d)  l^(cp) 

Here 

A,(^)  =     /  ^'  dq) 

J  J  / 


^      r  2    2 

■TX         vl"^'    cos  cp 


From  ( II-8),   we  get 


k'   cos  cp 


=     /    v^l-k  cos  ^    ^    K/l-k  cos  ^  j=±=: 

-rt  V  y  /i-k- cc 

=  y"        -/l-k'^cos^cp     1^     U/i-k'^cos^cp    |^Vj^,cp)J 


.^(v^+l)(k  sin''^  +  k'   sin  cp)V^(^,cp)  [^  ^  dcp 

\/l-k'   cos  cp 


dcp 
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In  the  last  equation,  (3-5)  'was  used  to  transform  the  integrand.   Next,  we 
integrate  the  first  term  of  the  integrand  twice  by  parts,  and  make  use  of  the 
2rt  periodicity  of  ^  and  V  .  The  entire  integral  on  the  right  hand  side  of 
(II-9)  becomes 


y  V^(^,cp)  <-  /l-k'^cos^cp  1^  f/l-k'^cos^cp  ^j 


-  V  (v  +1)  (k  sin  ■a  +  k'  sin  cp)  J. 
m^  m  '  ^  ^J 


dcp 


/    2   2 
\/l-k'  cos  q 


Using  (5-6),  this  becomes 


r                      r                           2   2n     ^1^^) 
y  V^(^,q))   [-n^-v^(v^+l)kSin^^J   ^^ 


/    2    2 
yl-k '  cos  q 


By  (II-8),  this  is  equal  to 


[^^.-v^(v^+l)k^sin^'fl]A  (-a) 


Thus,  we  have  shown  that  A.(^)  satisfies  (3-5): 


M^-o.%     |j  (v/l-i^'cos^,  ^_)  ^  [vJv^^Dk^sin^.^jA^C.)  =  0 
By  the  matching  condition  V  (0,rp)  =  -V  (0,-fp),  and  the  evenness  of 


^  (cp)  in  cp,  (II-8)  gives  us 


Aj(0)  =  0 


By  the  rigid  boundary  condition  V^^(^^,q))  =  0,  (II-8)  gives  us 


Aj(.^)  =  0 


Hence,  from  (II-IO),  (ll-ll),  (11-12) 


Aj(^)  =  Cjye(^;k,v^,Hj) 


2k 


Thus  every  term  A   {■&)   ^.(cp)  of  the  sum  In  (ll-?)  Is  an  even  Lame'' product 
In  (5.11). 

Finally,  there  can  be  only  a  finite  number  of  terms  in  the  sum  (II-7). 
First,  we  note  that  there  can  be  only  a  finite  number,  say  q,  of  independent 
terms,  since  each  product  in  the  sum  is  an  eigenfujiction  belonging  to  the 
eigenvalue  v  ,  and  an  eigenvalue  of  (5.5)  can  have  only  a  finite  multiplicity. 
Second,  there  can  be  no  more  than  q   terms  altogether,  since  any  q+1  terms  satisfy 
a  linear  relation 

q±l 

2^   CjAj(^)  J^    (cp)  =  0  (11-15) 

J»l 

But  we  can  find  some  i8  for  which  all  the  A.(i8)  are  different  from  zero;  for 
this  value  of  -Q ,    (II-15)  constitutes  a  linear  relationship  between  the  ^.(cp), 
which  is  impossible. 

If  V  (^,cp)  is  an  odd  function  of  cp,  the  proof  proceeds  in  the  same  way 
except  that  the  odd  Lame  fvinctions  are  used,  and  ( 11-11 )  is  replaced  by 
A  (0)  =  0  which  follows  from  the  other  matching  condition  V  (0,cp)  =  V  (0,-cp). 

The  above  proofs  of  the  orthogonality  and  completeness  of  the  Lame 
products  are  based  on  related  proofs  presented  in  [ll]  and  [lo]  repectively, 
for  the  elgenf unctions  of  a  complete  ellipsoid. 


[10]    pp.   1+12-411+ . 
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